Iranian Journal of Numerical Analysis and Optimization oz 


Vol. 9, No. 1, (2019), pp 69-92 5 
DOI:10.22067 /ijnao.v9il.68851 


The network 1-median problem with 
discrete demand weights and traveling 
times 


M. Abareshi and M. Zaferanieh* 


Abstract 


In this paper, the 1-median location problem on an undirected network 
with discrete random demand weights and traveling times is investigated. 
The objective function is to maximize the probability that the expected sum 
of weighted distances from the existing nodes to the selected median does not 
exceed a prespecified given threshold. An analytical algorithm is proposed to 
get the optimal solution in small-sized networks. Then, by using the central 
limit theorem, the problem is studied in large-sized networks and reduced 
to a nonlinear problem. The numerical examples are given to illustrate the 
efficiency of the proposed methods. 


Keywords: Facility location; 1-median problem; probabilistic weights; prob- 
abilistic traveling times. 


1 Introduction 


In the classical deterministic p-median problem, the purpose is to locate p 
new facilities on the links or nodes of a given network G, so that the total 
weighted distances from all nodes that are considered as demand points to 
the nearest facility, is minimized. The p-median problem on general net- 
works is NP-hard; see [8]. In a tree network, Kariv and Hakimi [12] proved 
that the problem can be solved in O(p?n?) time. Tamir [18] improved the 
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time complexity on tree networks to O(pn?). Hakimi [10] demonstrated that 
the classical 1-median problem has nodal optimal solution, and Goldman [9] 
presented a linear time algorithm to find the optimal solution. For the case 
p = 2, an O(nlogn) time algorithm was provided by Gavish and Sridhar 
in [7]. 

Determining precise demand weight of each node may be impossible. For 
example, the number of patients arising from different regions going into a 
hospital cannot be predicted precisely. Therefore, it is desirable to treat 
demand weights as random variables. Frank considered a probabilistic case, 
in which the demand weights are random variables with arbitrary known 
or estimated probability distributions, [5]. He discussed the median and 
center problems on a network with independent random demand weights 
and extended the classical 1-median problem to the maximum probability 
absolute median (MPAM) problem. The MPAM problem seeks to find 
a new facility maximizing the probability that the total weighted distances 
from all nodes to the new facility is less than or equal to a given threshold. 
He also defined a maximum probability center to maximize the probability 
that the maximum weighted distance is less than a prespecified threshold 
value. Later, the results for the maximum probability median were extended 
by Frank to the case where demand weights are correlated normal random 
variables, [6]. 

Berman and Wang [2] extended the problem to some single facility lo- 
cation problems as maximum probability 1-median, 1-anti-median, 1-center, 
and 1-anti-center problems with independent discrete random weights and 
presented efficient algorithms to solve these problems. In [3], They also stud- 
ied the p-facility location problem with maximum probability and formulated 
it as a linear integer programming problem, then determined the solutions in 
large-scaled networks by using the normal approximation. 

Puga and Tancrez [16] studied a location-inventory problem for the de- 
sign of large supply chain networks with uncertain demand and proposed a 
heuristic algorithm to solve the problem. Bieniek [4] considered the single 
source capacitated facility location problem with independent and identically 
distributed random demands. The unified of a priori solution for the loca- 
tions of facilities and for the allocation of customers to the operating facilities 
was found. 

The distances between nodes can also be considered as nondeterministic 
parameters. The uncertainty about the lengths of the links of a network often 
occurs when the lengths are measured in units of traveling times instead of 
geographical distances. Mirchandani and Odoni studied the p-median prob- 
lem when traveling times on the links are discrete random variables, [13]. 
The p-median problem with continuous random lengths of links, was studied 
by Handler and Mirchandani, [11]. They formulated the problem for locating 
p medians and provided an algorithm for the 1-median problem. The uncer- 
tainty in link traveling times can also be considered in some other location 
problems. Berman et al [1] studied the maximum covering location problem 
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on a network with traveling time uncertainty represented by different travel- 
ing time scenarios. The proposed models were applied to the analysis of the 
location of fire stations in Toronto, Canada. 


In this paper, a new case is proposed in which the weights of nodes and 
the traveling times on links are both discrete independent random variables. 
The problem is to find the location of the new facility in such a way the 
expected sum of weighted distances is less than a prespecified value with 
maximum probability. The problem is studied in both small and large-sized 
networks. In Section 2, the problem formulation is stated and an expected 
1-median with the maximum probability is defined. The solution approach 
for small-sized networks is stated in Section 3, and a precise algorithm to get 
the optimal solution is presented; then a numerical example is provided. In 
Section 4, by using the normal approximation, the solution is determined in 
large-sized networks and a numerical example is presented. Summary and 
conclusions are stated in the last section. 


2 Problem formulation 


Let G = (N,E) be a network, where N is the set of nodes and E is the set 
of links with |N| =n and |E| =m. The weight of each node v; EN, h= 
1,...,n is denoted by wy. Let the link connecting two nodes vu; and v; and 
its length be denoted by e;; and 1;;, respectively. In a deterministic network, 
the shortest distance between two points y,x« € G is denoted by d(y,x). Let 
x be a point on the link e;; at distance x from node v;; then the shortest 
distance between x and an arbitrary point y is given as follows: 


a(y,2) = min{d(y, v;) +2, d(y,v,;) +4; — 2}. (1) 


Hence the shortest distance from point y to x € e;; is covered via node 1; 
or node v;, whichever that is shorter. In a stochastic network, each link e;; 
has a deterministic geographical length |;;, that is a real positive fixed value. 
Moreover the traveling time along the link e;; is a random variable dependent 
on the traveling speed. In this case, the distances are computed in terms of 
traveling times. It is clear that the traveling speed along each link may be 
changed under such various conditions as traffic congestion. 


To introduce the traveling time of a link, we use the inverse of traveling 
speed. Let the traveling time and inverse of traveling speed along the link 
ei; be denoted by 7;; and SP,;;, respectively. If SP;; is a discrete random 
variable, T;; is also a discrete random variable equal to SP,j;l;;. In this 
paper, the random values for the inverse of traveling speed of the link e;; 
are considered to be independent from the other links and chosen from the 
set {spi;; SDi;, ashests spi;?} , where U;,; is the number of all states of traveling 


72 M. Abareshi and M. Zaferanieh 


speed in the link e;;. The relevant probability to spj; on the link e;; is Pj 
for u=1,...,Ui;. 

Because of random values for traveling times, the network G gets into 
different state-based networks G,. Let Gstate = {G,, r = 1,...,R}, where 


R= Newer U;;. We define 


WD,(a) = )° wad;(vn,2), r=1,...,R (2) 
h=1 


as the total weighted traveling time from all nodes to the point x in the rth 
state G, € Gstate, in which d,(vp,, x) is the distance between vp, and x € G,, 
that is determined by equation (1) noting that the traveling time of each 
link e;; € E is computed by spj,l;;. Then the expected 1-median problem is 
defined as below. 


Definition 1. The expected 1-median problem in network G with probabilistic 
distances seeks to find a point x which minimizes 


R 
E(WD(2)) = >> P.WD,(2), (3) 


where W D(x) is the random variable denotes the weighted traveling time from 
all nodes to the point x having the values WD,(x), r =1,...,R computed 
by expression (2), and P,. is the probability of the rth state in set Gstate. The 
function E(W D(a)) is known as the expected 1-median function. 


In continue, an independent discrete random weight W), is allocated for 
each node vu, € N. The discrete values of W;, are denoted by w,[k], k = 
1,...,A,, where K;, is the number of different situations for the weight of 
node vp, € N and it is assumed that w,[k] < walk +1] fork =1,...,K,—-1. 
The probability corresponding to wa[k] is denoted by Pi [k], k =1,..., Kn. 
Let T be a preselected admissible threshold. The proposed objective function 
is maximizing the probability that the expected 1-median function (3) does 
not exceed the given threshold 7. Here T is an upper level, which is given 
by the decision maker. 


Definition 2. Let the weights of nodes and the lengths of links be independent 
discrete random variables with known probabilities. A point x € G is an 
expected 1-median with the maximum probability if it maximizes the following 
probabilistic objective function: 


R 
PO P,WD,(e) <7) (4) 


and consequently 
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R n 
P'S" P. >” Wadr(vn, 2) <T), 
h=1 


r=1 
where P'(.) is the probability corresponding to the random weight vector of 


nodes. 


The expected 1-median with maximum probability is N P-hard on general 
networks; see [2]. 


2.1 Calculating a lower bound 


In this section, we provide a lower bound approximation for the objective 
function (4) at each point x of G. The lower bound can be calculated as an 
expected function, which would be readily calculated in comparison with the 
probabilistic objective function (4); see Lemma 1. 


Lemma 1. [f up = pas wn[k|P) [k], h =1,...,n is the expected value of 
W),, then the probability Pi, P,WD,(x) < T) is greater than or equal 
to 1— ee P, i Eh d,(Un, x). 
Proof. Based on Markov’s inequality, if X is a nonnegative random variable 
E(x), E(X) 
anda > 0, then P(X >a) < or in other words P(X < a) > 1-——., 
a 


a 
where F(X) is the expected value of X. Therefore, Since Wi,...,W, are 
independent random variables, the following inequalities hold: 


R n R n 
P'S" Pe >" Wide(on, 2) <T) > 1-EO OP.) del, x)) 
h=1 


r=1 r=1 h=1 
R n 
E(W; 

=1- OR a .2) 
r=1 h=1 
R n in 

=1- P,. —d,(vp, x), 
ZF Dap aew) 


which completes the proof. 


Lemma 1 determines a lower bound approximation for the function (4), 
when a given point is available as the solution for the expected 1-median with 
maximum probability problem. This helps to correct the available value of 
T and estimates a more reasonable amount to the preselected threshold. To 
determine the solution of the expected 1-median with maximum probability, 
we define the antipode points which are originally introduced by Frank, [6]. 
Using the antipode points, we bring forward the state-based primary regions 
in links, which are useful to estimate the distances; see [2]. 
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Definition 3. For each node vz, € N and each stater,r =1,...,R, the point 
a on the link e;; is called an antipode associated with the node Up in state 
r, if the following equation holds: 


2, Spi, + d,(vp, vi) = (lig — 2, )eDl, + d,(vn,v;)- 


An antipode z/,, on the link e;; provides two regions [v, 2/7 7) and [2 he ny Vil; 
where the fatanices between the points on the link e;; in each of these regions 
and the node va, in state r can be distinctly computed. Considering all 
antipode points and gathering them up together yield some basic intervals 
on e;;, which we call them the state-based primary regions; see Definition 4. 


Definition 4. The set of all points between two consecutive antipodes a 


and z ‘s in the link e4; Genceiaes with the node v1. € N in the state r’ and 
node uv; € N in the state r”, (r’,r” © {1,...,R}) ts called a state-based 


V, V, 
primary region Exe oh 


Recall that the state-based primary regions are special with a remarkable 
Properly that the distances from any node v;, in each state r to all points 
Ze (zi), zit, x] are computed either by xspj; + d,(vn,vi) or (lig — 2) spi; + 
d-(Un, ¥;) and it remains unchanged throughout this region. By using this 
property, the 1-median objective function in the state r, that is, WD,(a) = 
1 Wad-(vp, x), in each state-based primary region can be computed. 


Next, a reformulation of the expected 1-median problem with maximum 
probability is given. Since the weights of nodes are discrete random variables, 
there are different situations in network. Let the set of all situations for 
- vector of weights be denoted by Wstate = {W*|W®* = (wj,...,w3),s = 

., 5}, where S = J];_, K;, and W* denotes the vector of weights in 
ca s. Therefore, the weight of node v;, h = 1,...,n, in situation s, 
s = 1,...,5 and its probability are denoted by w7 and: ios respectively. 
For all W* © Wetate and « € G, the following characterization function is 
provided to calculate the objective function (4): 


1 yo Pe oh, wid, (up, 2) <T 
Yws = % r=1°T h=1 “her , ey 
ee { 0 otherwise. o 


By applying this characterization function, the objective function (4) 
would be rewritten as follows: 


max f (x >e Yws(2)Piys, (6) 


where Py. =[],, en P,*®; see [2]. 
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3 The solution approach 


In this section, two lemmas are provided to determine the optimal solution of 
the expected 1-median problem with maximum probability when the thresh- 
old T is selected from some regular intervals. Furthermore, when T is selected 
out of these regular intervals, an algorithm is represented to find the opti- 
mal solution; see [2]. The l-anti-median and 1-median problems are used in 
Lemmas 2 and 3 to indicate some regular intervals for threshold T, where 
all points in the network G are optimal with probabilistic objective function 
equal to 1 and zero, respectively. In 1-anti-median problem, the purpose is 
to find the location of a facility in the network, where the total weighted 
distances from all nodes to the facility is maximized. The point x is 1-anti- 
median, if it maximizes }7/_, wnd(vn,z), where the weights and distances 
are fixed values. For more details and the solution approaches, see [14]. 

Let x, be the optimal solution to the following 1-deterministic-anti- 
median problem: 


Se Wh [K7,]d>(vn, a) 


for r = 1,...,R, where w,[K;,] is the largest probabilistic weight corre- 
sponding to node vp. Also let x/. be the optimal solution to the following 
1-deterministic-median problem: 


n 


min S- wp[l]d-(un, 2), 


for r= 1,...,R, where w»[1] is the smallest probabilistic weight correspond- 
ing to node v_z. Lemma 2 provides an upper level value to the threshold T, 
while Lemma 3 provides a lower level value to it. If 7’ is selected greater 
than the estimated upper level or smaller than the estimated lower level, 
then the maximum probability expected 1-median solution would be readily 
determined. 


Lemma 2. Jf A = argmax;ai,...R >>,—1 Wn[Knld;(vn, tr) while x4 is the 
corresponding 1-deterministic-anti-median solution to the following problem: 


n 


Kyld 
ax) wal n|da(vn, 2) 


and T is given in such a way to satisfy the constraint 


>. wn[Kp|da(un, ta) < T, 
h=1 


then each point x € G is the expected 1-median solution with the maximum 
probability equal to 1. 
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Proof. The weights of the nodes vz, € N are assumed to be arranged in 
ascending order; therefore, the following inequalities hold: 


wll] < wh{2] ees wn [Kp]. 


Considering the above inequalities together with the definition of x, yields 
the following inequalities for an arbitrary «: 


3 


n 


n R 
< SOP. So wrlKaldr(vn,t) < S° Pp S> wr[Knldr(vn, 2r) 


1 h=1 r=1 h=1 


R n n 
< SOP. So wnlKalda(on,ta) = S> wn[Kalda(on,va) ST. 
r=1 h=1 h=1 


Hence Yy: (x) = 1 for all W* € Werate and x € G; therefore, the objective 
function value at every point x € G is equal to 1, that is, all x © G are 
optimal. 


Lemma 3. Let B = argmin;=1,....r >>,—1 Wallldr(vn, 2.), where x'y is the 
1-deterministic-median solution to the following problem: 


min >a |dp(un, 2). 


reGB 


fT < aa wall]da(un, 2’), then all points « € G are optimal to the 
expected 1-median problem with maximum probability equal to zero. 


Proof. Considering the definitions of x’,, the following inequalities hold: 


r <ul [1]dp( (vpn; 2p) se Ya [1]dp( (vn, 2p) 


h=1 
R n 
< So P. So wallldr(on, x N<OP, Sal r(Up, 2) 
r=1 
R 
<K > P,S wnylk]d,(up, 2). 


The foregoing inequalities yield Yy«(a) = 0 for all a € G and W® © Wetate; 
therefore, all points « € G are optimal with probability equal to zero. 


Lemma 4 provides a sufficient condition to ignore some links, which along 
them, the characterization function (5) is equal to zero. 
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Lemma 4. Let v; and v; be two adjacent nodes, where Yws (ui) = Yws (vj) = 
0, for some W* € Wetate; then Yws(x) is also equal to zero for all x € e;;. 


Proof. By the assumption, Yws(v;) = Yws(v;) = 0; therefore, the following 
inequalities hold: 


R n 
ys P,. a w7,dr(Un, v;) > T. 


Furthermore, each function d,(vp,2) is concave in each link e;;; see Figure 
1. Therefore, the function: 


R 
S- P, ye w;,d,(Up, 2) 
r=1 h=1 
is concave on the link e;;; see [14]. This yields the inequality 


R n 


S> P, x wrdr(un, 2) > T 


r=1 h=1 


for all x € e;; and consequently Yw:(x) = 0, which completes the proof. 


dy (UpsX) » 


de(p, 04) 


d, (Up, 0) 


0 ti spijlij 
i 
Zen 


Figure 1: The figure of function d;(vp, x) on the link e;;. 


Let the threshold T be selected between the upper and lower bound values 
provided by Lemmas 2 and 3, 


S> wall]da(vn, 2p) <T< S> wn[Knlda(vn, 2a). 
h=1 h=1 
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In this case, the solution of the expected 1-median problem with maximum 
probability would be determined by defining some new points x in network 
links, where the characterization function Yyws(#) for the weight vectors 
W®* © Wetate would change. First some special points in stated-based pri- 
mary regions, called jump-points, are specified as introduced in [2]. Then it is 
shown the objective function (6) would be changed only at the jump-points. 
The shortest path from the points in a state-based primary region on the 
link e;; to each node vp, would be obtained through the node v; or v; without 
switching to the other node. Hence, the value of objective function (6) can 
be readily computed in all jump-points in each state-based primary region. 
Therefore, the optimal solution in each primary region would be determined. 
By comparing all solutions in the state-based primary regions together, the 
local optimal solution in each link would be obtained. 


Let [a wD z, ,,| be a state-based primary region associated with the nodes 
and & in the link e;; in states r’ and r”, respectively. Then the set of nodes 
N and the set of states corresponding to the traveling times of links, G'state, 
can be partitioned into two sets L, and K, ay as follows: 

La fire, = a ahs Ken = {(h, r)lze, = a gh 
Indeed, the set L,, includes the pair indices (h,r), where h corresponds to 
the node va € N and r corresponds to G; € Gstate, so that the position 
of antipode es in the link e,;; is before the position of the antipode re 
associated with the node wv in the state r’. The set oa is defined similarly. 
Therefore, for (h,r) € L,,, the distance between any point x € ees er) 
and the node vp, in the state r is determined through the node vj, that. is, 
dy-(Un, x) = spi; (liz — 2) +d,(vn, vj), where sp, is the inverse of the traveling 
speed on the link e;; in the state r. The corresponding value for d,(up, 2), 
where (h,r) € K+, would be determined through node v; in a similar manner. 


As an example, consider the small network shown in Figure 2 with three 
state-based networks G,,G2, and G3 with spi, = 0.5, spi, = 0.3, spo3 = 
0.3;. sha = 0.6 tor P= 19.3 and. spi, = ep, = 0:2, epi = apy = 03 
and sp}, = sp3, = 0.4. Consider the link ej2, then the set of antipodes 
associated with the node 3 is {2}%, 293, 733} = {3.7, 2.7, 2.7} and with the 
node 4 is {21%, 23%, 23%} = {3.4,4.1,4.8}; see Figure 2. For the primary 
region (23%, 21%] = [2.7,3.4], we have Ly = {(2,1), (2, 2), (2,3), (3, 2), (3,3)} 
and Ki = {(1,1), (1,2), (1,3), (3,1) (4,1), (4,2), (4,3)}. Note that for all 
states r, the antipodes corresponding to the nodes 1 and 2, are the points 
lig = 6 and 0, respectively. Now consider one of the weight vectors W* € 
Wstate with components w,h = 1,...,n. So for any point x € [277 ), 271g], 
the following relationships hold [2]: 


S > Pe >) widr(vn, 2) <T 
r h 
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Figure 2: The state-based primary regions in the link e2 in a small-sized network. 


=> S- P,w7,d,(Upn, £) + S- P,w7,d,(Up, £) < T 


(r,h)eELy, (r, nek, 
or 
s Pw (dr (vn, vj) + spj; (lig =a) + ‘> Pw (dr(vn, vi) + Spj;2) ay i 
(r,h)eLy, (r,h)EKs, 
=> S- Pw}, (dr(vn, v;) + sp; liz) + S- P,wpdr (Un, Vi) 
(r,A)EL., (r, nek, 
Soa S- P,w),8Pi; — S- P,w7,spjj)t < T. 
(hr)eKs, (h,r)ELy, 


(7) 


Definition 5. If inequality (7) is binding at a point «™°, then x” is called 
a jump-point with respect to weight vector W* © Woetate- 
If we set 
S- P,wy, (dr (Un, Vj) + SP} jli3) = AW: > P,wydr (Un, vi) = BY 

(r,h)ELy, (r,AEeKS, 

and 
we 

ee er ee ee 
(h,r)EK*, (hr)eL, 
then ee we 
Pia _ ‘Ee — AE = Bri Kk 


Ws 
COL), (0K) 
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Let J be the set of all jump-points in state-based primary region [2d wD Bor: tl 
unified with two boundary points {2/7 ,, 2/7, ;}- 

Lemma 5. The expected 1-median objective function would change only in 
the points belonging to the set J. 


Proof. The objective function value would change in some point z, if the 
value of Yys(x) is changes from 0 to 1 or vice versa for some vectors W*. 
This guarantees that these points belong to the set J. 


If a local optimal solution in a state-based primary region in some links is 
given, Lemma 6 provides some optimal intervals in this region. The coefficient 
of variable x in inequality (7) is a criterion, which determines the optimal 
intervals; for more details see [2]. 


Lemma 6. Let «”° € ba #7 | be the local optimal jump-point while yW 
and yw? are the left and right consecutive jump-points to a", respectively. 
Then the optimal interval in [ee is zen pl is determined as follows: 
s 1 1 s 
1. If Cle by (rt) > 0 and Ce 5) (rt) <0, then the interval [y ,x”"] 
is optimal. 


2. If Ort) (or) <0 and OF ici > 0, then the interval [xW*  yV*) 
is optimal. 


&. If Cle BE) e!",K) >0 anda” is the first jump-point after aes then the 
interval (21 ,,a”"] is optimal, where the starting point of the interval 
can be distinctly verified. 


4. If Chr) (r’,K) <0 and «™” is the last jump-point before es then 


the interval ee, x) 18 optimal, where the end point of the interval 
can be distinctly verified. 


Proof. The coefficient of variable x in inequality (7) is OM ny cet, x) lf 

s ue s 
Ch by (9K) 2 0 and Cle 1) (tt Ke) = 0, then Yws (x) = 1 forall x < ev 
and Yywi(z) = 1 for all « > y[” . For other vectors W’, which their 
respective jump-points x lie outside the interval [y!”’,r“°], the values of 
Yw-(a) remain unchanged in this interval. Note that if OP pyc > 0 and 
Ce ry (eK) > 0, then Yyi(2"") = 0. Therefore, f(z¥") < f(y"), which 
contradicts the optimality of point «”°. The other cases can be proved 
similarly. Note that the rest of states for COM 1) (r",K) Coe ny tnt K)? and 


On a (r”,K) are impossible, since they contradict the optimality of point 
ws 


x 
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3.1 The Algorithm: finding a global optimal solution 


In this section, Algorithm 1 is represented to describe that how a local opti- 
mal solution can be obtained on an arbitrary link and consequently a global 
optimal solution is determined. Recall that if two arbitrary adjacent nodes, 
say vu; and v;, have the same zero objective function values, then all points on 
the link e;; have also objective values equal to zero. Therefore, these edges 
should be ignored in the proposed algorithm. In other side, as described in 
Lemmas 2 and 3, the solution can be determined directly for spacial values of 
T. The idea of Algorithm 1 is to specify the state-based primary regions on 
all links and obtain the local optimal solution in all regions. Hence, G'state the 
set of all state-based situations for the network G, , and the set of all weight 
situations, Wstate, are determined. Then the primary regions are located 
by computing the antipodes while the set of jump-points are determined in 
each primary region, as well. The corresponding value of the objective func- 
tion (6) is estimated for all jump-points x by evaluating Yw:(a) for vectors 
W* © Wetate. By comparing these solutions, the optimal solution in each 
primary region and consequently on each link would be determined. Finally, 
the global optimal solution in the network would be determined by compar- 
ing the local optimal solutions on all links, while the optimal regions are fully 
specified by Lemma 6. 

In the represented algorithm, the maximum number of state-based pri- 
mary regions in each link is (n — 2)R + 1 and the maximum number of 
jump-points in each link is S((n—2)R+1), where S and R are the number of 
states for the weights of nodes and the traveling times of edges, respectively. 
To specify the solution, the values of Yws, s=1,...,S for each jump-point 
should be estimated. Therefore, if any link is not ignored by Lemma 4, in the 
worst case, the time of algorithm in network would be mS$?((n — 2)R +1). 
This time would be so long, particularly for large-sized networks. In section 
5, we propose a method to solve the problem in large-sized networks without 
using the Algorithm 1. 


3.2 Numerical example 


In the following, a numerical example is provided. The optimal solutions are 
estimated for different values of threshold T by performing the steps of Algo- 
rithm 1 in MATLAB R2014a, and the sensitivity analysis for different values 
of weights is proposed. Using Lemmas 2 and 3, some regular intervals for 
threshold T are introduced, where the optimal solutions along them would 
be readily obtained. 


Example 1. In this example, we study an squared network given in Figure 
3 with 25 nodes and 40 links. The lengths of the links have been written next 
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to the links while the weights of the nodes with their relevant probabilities 
are shown in Table 1. The traveling speeds in different states for each link 
ej; with their relevant probabilities are given in Table 2. 


Algorithm 1 Determining the expected 1-median with maximum probabil- 
ity in the link e;; 


1: Determine the solution of the deterministic prob- 
lems min,=1,...,.R Ming eG, >-,—1 Wall|dr(vn, 2’) and 
MaX;=1,...RMAXzeG, >, —1 Wn[KnJd-(vn,z), set their optimal val- 
ues as I and U, respectively. If T < L or T > U, then the optimal 
solutions are determined by Lemmas 3 and 2, respectively. Otherwise go 
to step (2). 

2: Consider the set of all state-based situations for the network G' as Ggtate = 
{G,, r =1,...,R} in which R = Ilex U;; with relevant probabilities 
P, = Il.j,en Pj, where Pj; is the probability that the link e;; has the 
traveling time d,(v;,v;) = Spi jliz. 

3: Arrange the demand weight situations in the set Wstare = {W*| W* = 
(wj,.--,W,),8 =1,...,S} with relevant probabilities Phys = [[,, cn Pe 
for s=1,...,8. 7 

4; Determine and sort the set of ZY = {z,|h=1,....n, h#ij, r= 
1,...,£t} U {0, f,,} on the link e,; € # in ascending order. Rename the 
members of Z” as zy,k =1,...,|Z”|. 

5: Set Wate = state\{W*|Yws (v;) = Yws (v;) = 0} and k=1. 

6: Determine the state-based primary regions [z,, z,.41], where zp, 241 € 
Zi, and the jump-point set J, = faklak = gM ec 
[2k Zk-+1| for some vector W* € Woes d= 1: .,D}, where D = | Jz|. 
Set. = 1, = Land f(e8) = 0. 

7: If Ywa(ak) = 1, where W% is the qth member of W%,,;., set f(a‘) = 
f(a*) + Piya, else go to step (8). 

8: Setq=q+l1. Ifq< |W.44-|, then go to step (7), else set d= d+ 1 and 
go to step (9). 

9: If d< D, set f(x*) = 0, ¢ = 1, and go to step (7), else go to step (10). 

10: Set vk, = arg maxg f(x). 
11: Setk=k+1, ifk <|Z¥|—1, go to step (6), else go to step (12). 
12: Set re,, = aTGMax_=1 (aed ee) 


pivwy 


The upper and lower bound values provided by Lemmas 2 and 3 are 
5716.435 and 310.2833, respectively. Therefore, if T > 5716.435 or T < 
310.2833, then all points of network G are optimal with probabilities equal 
to one and zero, respectively. In addition, using the proposed algorithm, it 
can be numerically verified that, for values of T’ < 1050 and T > 1900, the 
optimal objective function value is equal to zero and one, respectively. 

To verify Algorithm 1 in details, we consider the threshold T = 1100. By 
applying Lemma 4, only the links €7,12, €8,13, €11,12, €12,13, €12,17, €13,14 and 
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Figure 3: A squared network with 25 nodes. 


Table 1: The weights of nodes in different states 


nodes I 2 3 4 5 6 7 8 9 10 11 12 13 
state 1 75 35 33 59 47 57 56 72 23 86 78 62 53 
state 2 69 68 18 88 74 54 28 5 90 90 40 50 96 
state 3 44 36 4 67 9 50 57 = 25 72 =i 32 55 4 
state 4 47 94 57 28 80 43 82 27 7 38 50 78 39 
state 5 75 28 33 76 43 33. 5372 60 44 41 15 68 
nodes 14 15 16 17 18 19 20 21 22 23 24 25 Probability 
state 1 82 18 51 75 48 27 «672 8 32 46 74 4 0.2 
state 2 77 18 10 +54 6 9 25 88 51 1 66 73 0.25 
state 3 57 6 19 16 69 23 41 27 8656 51 63 62 0.15 
state 4 101 32 38 97 107 46 90 68 7 28 24 87 0.1 
state 5 47 17 36 9 57 29. 58 44 48 66 72 48 0.3 


Table 2: The different traveling speeds 


State 1 2 3 4 5 6 


link traveling speed | 100+i+ 7  7O+i+ 7  50+i+7 40Fi+7 S8OFi+7 90FI4j 


probability 0.15 0.2 0.1 0.25 0.23 0.07 
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€13,18 have to be searched for state s = 3 and the other states s = 1,2,4,5 as 
well as the other links would be ignored. We follow the steps of algorithm for 
links €j213 and €13,14. The antipodes of the link e;213 are sorted ascending 
according to their distance from vertex 12, as follows: 


Z'*13 — £0, 5.5061, 5.6330, 5.8108, 5.8753, 5.9289, 5.9742, 9}. 


The jump-points obtained by equality (7) are all outside the primary regions 
and the only jump-points of the link e12 13 would be its antipodes. Computing 
the values of Suan P,. Sar w7d-(Un, Zz) for z, € Z'*13, k =1,...,8 results 
in Yws(z~) = 1. So all points in the link e12,13 are optimal with the objective 
function value equals to 0.15. The other evaluated link is e13,14, where its 
antipode set is 


71314 — (9, 1.5085, 1.6971, 1.9164, 2.1749, 2.8602, 3.3285, 3.8673, 5.3977, 
7.4310, 8.1409, 8.7211, 9.2043, 33.5780, 34.0291, 34.6129, 
34.8126, 34.9741, 35.1074, 44.9147, 45.7261, 45.7886, 46.9351, 
46.9441, 47.4125, 47.8139, 48.0984, 48.1616, 48.2600, 48.5373, 
48.6747, 48.9968, 49.1737, 49.2538, 49.4117, 
49.6128, 49.7849, 69.6898, 71.4483, 72} 


Again, the jump-points obtained by equality (7) are outside the primary 
regions and the set of jump-points includes just the antipodes. For antipodes 


{0, 1.5085, 1.6971, 1.9164, 2.1749, 2.8602, 3.3285, 3.8673, 5.3977} c Zi314 


the value of Yys is equal to 1. Therefore, the optimal interval obtained by 
Lemma 6 is [0 , 5.3977] with optimal probability 0.15. Using Algorithm 1, the 
solutions for different values of T with their objective (Obj) function values 
are determined and inserted in Table 3. The obtained optimal intervals (Int) 
in corresponding links with running times (t) are also reported. Moreover, 
the variations of the objective function (4) with respect to different values of 
T are released in Figure 4. 


Next, we develop a sensitivity analysis to the solutions by certain pertur- 
bations in the weights of nodes and traveling speeds of links. The weights of 
all nodes are assumed to be equal to 75 in all 5 states with relevant probabil- 
ities {0.55, 0.1, 0.1, 0.1, 0.15} while the lengths of all links are considered 
equal to 50. In addition, 6 states are considered for links’ traveling speeds that 
are assumed to be the same for all links equal to {100, 70, 50, 40, 80, 90} 
in states r= 1,...,6 with relevant probabilities given in Table 2. The weight 
of node 1 is multiplied by 100 in the first state, that is, wt = 7500. This 
alteration results in the solution of maximum probability expected 1-median 
problem for T = 3700 being attained in node 13 with probability 0.45. More- 
over, Yws(a*) = 1 for the states s = 2,3,4,5 and Yyw:(a*) = 0, where a* 
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denotes the optimal solution, also the interval [0, 6.8090] in the link e13,18 is 
optimal. 


By increasing the value of T from T = 3700 to T = 6000, the optimal 
value of the objective function remains unchanged while the optimal solutions 
would be spread in the network. For T = 6000, all nodes except 1, 5, 21 
and 25 lie in the optimal solution set. If 7 = 6100, only the node 1 is 
optimal with objective function value 1. In this case, the value of T is larger 
than the expected sum of weighted distances from all nodes to node 1 in all 
states s = 1,2,3,4,5. Notice that if the median is in node 1, the expected 
sum of weighted distances from all nodes to node 1 is the same in all states 
Sibi: 


Next, we multiply the weight of node 25 in state 5 by 100 while the other 
information remain unchanged. For T = 3700,...,6000, the optimal value 
of the objective function is 0.3. The solution for T’ = 3700 is node 13, by 
increasing the value of T to TJ’ = 6000 the value of optimal probability does 
not change but the number of optimal points is increased. For T = 6000 all 
nodes expect 1, 5, 21 and 25 lie in the optimal points’ set. For T = 6100 
the optimal probability is changed to 0.85 whereas the optimal solution set 
is reduced to node 1. 


If the probability of the first and fifth states are changed together, that 
is, Py = 0.15 and Pj; = 0.55, then the solution would be in the last node 
with optimal probability 0.85. Also, if the first and the fifth states of nodes’ 
weights have just the same probability equal to 0.35, then for T = 6100, both 
nodes 1 and 25 would be optimal with optimal probability 0.65. 


+ Optimal objective function |: 
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Figure 4: Variations of objective function (4) for different values of T 
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Table 3: The solutions for different values of T 


T = 1100 link €12,13 €13,14 t = 0.51 
Obj =0.15 | Int [0 , 9] [0 , 5.3977] 
T = 1150 link €8,13 €11,12 €12,13 €12,17 €13,14 t = 0.49 
Obj =0.15 | Int | (28.8133, 41] [34.9193 , 42 (0, 9] (0, 11.0125] [0, 9.2043] 
T = 1200 link €8,13 €11,12 €12,13 €12,17 €13,14 €13,18 t = 0.66 
Obj =0.15 | Int | [8.5634, 41] [32.5787 , 42 (0, 9] (0, 13.3575] (0, 9.2043] (0 , 71] 
T = 1250 link €2,7 €6,7 €7,8 €7,12 €8,9 €8,13 €11,12 e12,13 | t=1.51 
Obj =0.15 | Int | [6.8588 , 9] [65.7371 , 68 (0 , 3] (64.0216 77] [0 , 4.1260] [0 , 41] [0 , 4] (0, 9] 
T = 1250 link €12,17 €13,14 €13,18 €14,19 €17,18 €18,19 €19,20 €20,25 | t= 1.51 
Obj =0.15 | Int , 13.3575] [0,, 9.2043] [0 , 30.9657] [0 , 33] (71.4691 , 76] [0 , 11.1206] [0 3.3058] 0 
71.4483 , 72 (66, 8237 , 70] 
T = 1300 link €12,13 €13,14 t = 2.29 
Obj =0.65 | Int | [7.0710 , 9] [0 , 2.1749] 
T = 1350 link €7,12 €11,12 €12,13 €13,14 t= 3.25 
Obj = 0.65 | Int | [74.6709 , 77] [38.7520 , 42 [0 , 9] [0 , 7.4310] 
T = 1400 link €13,14 t = 3.66 
Obj =0.9 | Int 0 
T = 1450 link €12,13 €13,14 t = 4.58 
Obj =1 | Int | [5.5061 , 9] [0 , 2.8602] 
T = 1500 link €12,13 €13,14 t = 5.66 
Obj =1 | Int [0 , 9] [0 , 9.2043] 
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4 Central limit theorem 


In this section, the expected 1-median problem with maximum probability in 
the large-sized networks is considered. Considering all state-based primary 
regions, determining the jump-points, and finding the value of objective func- 
tion in all of them make a large number of computations to be done that 
practically is useless. Therefore, we state an approximation method, which 
provides a near optimal solution to the expected 1-median with maximum 
probability, the approach proposed in [3]. 

Let the mean and variance of W),, which is the independent discrete ran- 
dom weight of node vw, € N, be denoted by pz, and o7, respectively. Using 
Lemma 7, the mean and variance of term sit P,. Se eee Wrd-(vn, £) 
are determined. Then by applying the central limit theorem, the problem is 
reduced to a nonlinear fractional programming problem. 


Lemma 7. Based on the central limit theorem, the term 
pet. Ree het...n Wadr(vn, ©) in large-sized networks has the normal 
distribution with the expected mean value 


Hewo= >, Pr So pndr(vn,2) 
h=l1sn 


TH1 pak 
and variance ow p = Dpaty..n P2 Cretan MEOH: 2) 


Paina P, weet cd Widr(un, x) — HEWD 


SEWD 


Therefore, ~ N(0,1), which 


yields the following: 
P’ So OP SS Wad-(un,2) <T) 
r=1,..58 h=1,...,n 


Die ce P, iW Wid,(vn, x) — KEWD Z T— ee) 
OEWD a OERWD 


=P 


T — WeEwp 

CEwD 
where ¢ is the cumulative distribution function of the standard normal dis- 
tribution. By considering the ascending property of ¢, the main problem can 
be transformed into 


i SB P, a ae Lndr (Un, ©) 
max (8) 


7 yf otek i i ae 07,42 (Un, x) 


The above nonlinear fractional programming problem can be solved by 
using the existed approaches such as parametric or direct methods; see [17]. 
But here we use a method that considers the problem in three cases. Let x* 
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be the optimal solution to the problem (8); then three cases can be occurred 
associated with T — pwp(x*) that are discussed in details below; see [3]. 


1 
1. First case: T — pewp(2*) > 0. Setting T — wewp(x) = -, the 
Y 


problem is transformed into the following problem: 


min DY) Pr DD onr?de(vn, 2) (9) 


TH1 pg R BEL, 21.97% 


st. Ty - S- P, S- Lnydy(up, £) = 1. 
r=1 R h=1 n 


ga 


To solve the above problem, the optimal point x in each state-based 
primary region should be determined. Let [z/7 ,, z/7, ,] be a state-based 
primary region associated with the nodes | and k in states r’ and r” 
in the link e,;, respectively. Then for all x € [a4 ‘3 zg, ,| and an ar- 
bitrary state r and node vp, either d,(vupn,2) = d,.(vp, Vi) + spij,x or 
dy-(Up,&) = dy (vp, vj) + spi; (liz — x). Substituting d,(vn, x) in problem 
(9) and setting ya = , the problem can be rewritten as a second order 
problem with one linear constraint in each primary region, which can 
be solved by using such nonlinear convex programming algorithms as 
Lagrangian dual method or linearization approaches; see [15]. Compar- 
ing the objective function values of the optimal solutions in all intervals 
yields the global optimal solution in each link and consequently in the 
network. 


2. Second case: T — uzwp(«*) = 0. The second case is considered 
when the first one is infeasible. In this case, it is sufficient to solve the 
following problem: 


min ) 
zx 
ea be 


Pr >) pnde(un,2). (10) 
Goengdt. h=1 n 


Using the concavity of the objective function as mentioned in the proof 
of Lemma 4, this problem has at least a nodal optimal solution. So it is 
sufficient to search among the nodes of the network. If the optimal value 
of the objective function (10) equals to T, then the obtained solution 
is optimal to the fractional programming problem (8) with the optimal 
value equal to zero. Otherwise, the third case should be considered. 


3. Third case: T—jsgwp(«*) < 0. If the first and second cases are infea- 
sible, this case is considered similar to the first one, but the constraint 
is changed as follows: 
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> P,. S- Lnydy(up,, 2) — Ty = 1, 


and the objective function (9) should be maximized. The local op- 
timal solutions of the resulted problem can also be obtained in each 
primary region by using such nonlinear programming algorithms as the 
Lagrangian dual method or linearization approaches; see [15]. 


Example 2. To illustrate the efficiency of the proposed method based on cen- 
tral limit theorem, a test problem from the address (http : //people.brunel. 
ac.uk/mastjjb/jeb/info.html) with 100 nodes and 396 links is considered. 
Four scenarios with probabilities {0.1, 0.25, 0.45,0.2} for the network links’ 
traveling speeds have been considered. For each link e;; the traveling speed 
in state r = 1, 2,3,4 is assumed as sp}, = (r +i+ 7). 


We have considered 1000 different situations for vector of weights with 
discrete uniform distribution and each of components for random weights of 
nodes are taken from the interval [0,850]. Using the calculated mean and 
variance for W,,h = 1,...,100, the proposed method based on central limit 
theorem is applied for different values of threshold T = 10+(1300 + 50i);i = 
1,...,10. In each case, the resulted nonlinear problems are solved by using 
the interior-point solver in Matlab R2014a, and the best optimal solutions are 
determined. Then, the optimality of the solutions is verified by comparing 
e — Hewp(2*) 

Hw D(2x*) 
limit theorem with the values obtained from the objective function (6) for 
the obtained solutions x*; see Figure 5. 


the values of estimated objective function @ ) by the central 


In addition, the variations of the estimated relative error 
T — pew d(2*) 
ae ee al 
THwD(2*) 

f(x*) 
solution x* and T is selected from the set {104(1300 + 507);7 = 1,...,10} 
is depicted in Figure 6. As it is seen, increasing the value of T results in 
increasing the objective function optimal value and therefore decreasing the 
relative error. 


o( 
, where f(2*) is computed by (6) for the optimal 


90 


; T 
Figure 5: Comparison between ¢( 


Estimated relative error 
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Figure 6: Variations of relative error for different values of T. 


5 Summary and conclusions 


The median problem is one of the most common problems in location the- 
ory specially with nondeterministic parameters. In this paper, the 1-median 
problem on a network with independent discrete demand weights and trav- 
eling times is investigated. The objective function is devoted to maximizing 
the probability that the expected 1-median function does not exceed a given 
threshold. First, a precise algorithm to obtain the optimal solution in small- 
sized networks is presented. Next, by using the central limit theorem, an 
approach to find the optimal solution in large-sized networks is proposed, 
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where the original problem is reduced to some linear and nonlinear problems 
with a linear constraint. The numerical examples are given to illustrate the 
efficiency of the proposed methods. 
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